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Abstract: We study confinement-deconfinement phase transition for heavy quarks in a soft wall holo¬ 
graphic QCD model. We consider a black hole background in an Einstein-Maxwell-scalar system and 
add probe open strings to the background. Combining the various configurations of the open strings and 
the phase structure of the black hole background itself, we obtain the confinement-deconfinement phase 
diagram for heavy quarks in the holographic QCD model. 



Contents 


1 Introduction 1 

2 Einstein-Maxwell-Scalar System 3 

2.1 Background 3 

2.2 Phase Structure of the Background 5 

3 Open Strings in the Background 8 

3.1 Wilson Loop 9 

3.2 Configurations of Open Strings 10 

3.3 Heavy Quark Potential 13 

4 Phase Diagram 16 

4.1 Confinement-deconfinement Phase Diagram 17 

4.2 Meson Melting in Hot Plasma 18 

5 Conclusion 19 


1 Introduction 

Confinement-deconfinement phase transition is an important and challenging problem in QCD. Near the 
phase transition region, the interaction becomes very strong so that the conventional perturbation method 
of QFT does not work. For a long time, lattice QCD has been the only method to study strong interacted 
QCD. Although lattice QCD works well for zero density, it encounters the sign problem when considering 
finite quark density. See [1, 2] for a review of the current status of lattice QCD. Recently, using the 
idea of AdS/CFT duality from string theory, one is able to study QCD in the strongly coupled region by 
studying its weakly coupled dual gravitational theory, i.e. holographic QCD [3-19]. In [20], we considered 
a Einstein-Maxwell-scalar system and studied its holographic dual QCD model. We obtained a family of 
analytic black hole solutions by the potential reconstruction method. By studying the thermodynamics of 
the black hole backgrounds, we found a phase transition between two black holes with different size. We 
interpreted this black hole to black hole phase transition as the confinement-deconhnement phase transition 
of heavy quarks in the dual holographic QCD model. 

On the other hand, the heavy quark potential is an important observable relevant to confinement. It 
has been measured in great detail in lattice simulations [21] and the results remarkably agree with the 
Cornell potential [22] 

V (r) =--+ asT + C, (1.1) 

r 
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which is dominant by Coulomb potential at short distances and by linear potential at large distances with 
the coefficient as defined as string tension. In QCD, the heavy quark potential can be read off from the 
expectation value of the Wilson loop along a time-like closed path C, 

(W(C)) ~ (1.2) 

In string/gauge duality, the expectation value of the Wilson loop is given by [23] 

{W {C)) = J DXe-^^^, (1.3) 

where Sng is the string world-sheet action bounded by the loop C at the boundary of an AdS space. In 
[24-35], a probe open string in an AdS background was considered. The two ends of the open string are 
attached to the boundary of AdS background and behave as a quark-antiquark pair. Thus the open string 
could be interpreted as a bound state, i.e. meson state, in QCD. By studying the dynamics of the open 
string, the expectation value of the Wilson loop can be obtained, so as the heavy quark potential. From the 
behavior of the heavy quark potential, one is able to study the process that an open string breaks to two 
open strings with their two ends attaching to the AdS boundary and the black hole horizon, respectively. 
This string breaking phenomenon describes how a meson melts to a pair of free quark and antiquark in its 
dual QCD. 

In this work, we put probe open strings in the background obtained in [20]. We study the dynamics of 
the open strings to obtain the expectation value of the Wilson loop as well as the heavy quark potential. 
In [20], various black hole phases for different temperatures have been obtained. In this work, we found 
three open string configurations for the various black hole phases as in figure 1. According to AdS/QCD 
duality, these different open string configurations correspond to the confinement and deconfinement phases 
in QCD, respectively. This supports our preferred interpretation that the black hole to black hole phase 
transition in the bulk corresponds to the confinement-deconfinement phase transition of heavy quarks in 
the dual holographic QCD in [20]. Nevertheless, we found that the phase transition temperatures obtained 
from the black hole phases and the open string configurations are not exact the same. In fact, we will 
argue that neither the black hole phases nor the string configurations alone could explain the full phase 
structure of the confinement-deconfinement phase transition in QCD. The string configurations tells us 
that whether the system is in confinement or deconfinement phase, while the black hole phase transition 
determines the location of the phase boundary. By combining the two effects together in this paper, we 
find a more natural picture to describe the phase diagram of the confinement-deconfinement transition for 
the heavy quarks in QCD. Furthermore, in the deconfinement phase, we also study the meson melting 
process by studying the process of an open string breaking to two open strings. 

The paper is organized as follows. In section II, we consider an Einstein-Maxwell-scalar system. We 
review how to get the analytic solutions in [20] by potential reconstruction method and study the phase 
structure in the these backgrounds. In section III, we add probe open strings in our black hole background 
to study their various configurations. We calculate the expectation value of the Wilson loop and study the 
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(a) no BH 


(b) small BH 


(c) large BH 


Figure 1. Three configurations for open strings in a black hole. There is no black hole in case (a), open strings 
are always connected with their two ends on the AdS boundary. In the small black hole case (b), open strings can 
not exceed a certain distance from the boundary and are still connected with their two ends on the AdS boundary. 
In the large black hole case (c), the open strings with their two ends far away enough will break to two open strings 
with their two ends attaching the AdS boundary and the black hole horizon, respectively. 


heavy quark potential. In section IV, by combining the background phase structure and the open string 
breaking effect, we obtain the phase diagram for confinement-deconfinement transition. We further study 
the meson melting process in the deconfinement phase. We conclude our result in section IV. 


2 Einstein-Maxwell-Scalar System 

In this section, we review the black hole solution and its phase structure obtained in [20]. 

2.1 Background 

We consider a 5-dimensional Einstein-Maxwell-scalar system with probe matters. The action of the system 
has two parts, the background part and the matter part, 


S 


Sb + S„ 


( 2 . 1 ) 


In Einstein frame, the background action includes a gravity field a Maxwell field and a neutral 
scalar field (j), while the matter action includes a massless gauge fields A^, which we will treat as probe, 
describing the degrees of freedom of vector mesons on the 4d boundary. 


Sb = 


IGttGs 


Sm=- 


J (fxy/^ 




( 2 . 2 ) 

(2.3) 


where is the coupling constant for the gauge field strength dfj,Ay — dyAfj_, f (cp) is the gauge kinetic 

function associated to the Maxwell field and V {(p) is the potential of the scalar field cp. 
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The equations of motion can be derived from the above action as 


^ (f2 + F^) , [f{ct>)Fn = 0, [/(</-)T7] = 0, (2.4) 




F^lpRv + 2 


d^4>di,4> - -g^y {d4>f - g^yV 


(2.5) 


To solve the background of the Einstein-Maxwell-scalar system, we first turn off the probe gauge field 
and consider the ansatz for the metric, scalar field and Maxwell field as, 


j 2 e 

as = — 


2A{z) 


2 dz^ 


—g{z)dt'‘ + + dx^ 

9{z) 

4> = (j){z), A^ = At {z ), 

which leads to the following equations of motion for the background fields. 


+[9_ + 2>A' + 


z'^e 


2g 


z^g 


= 0 , 


A'! + 


!L + a'-^]a', = o, 


k'2 




A" - + -a' + — = 0, 

z 6 

- e-^^z^fA!^ = 0 , 




= 0 . 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


To solve the above equations of motion, we need to specify the following boundary and physical conditions: 


1. Near the boundary z ^ 0, we require the metric in string frame to be asymptotic to AdS^; 

2. Near the horizon z = zh, we put the regular condition At {zh) = g (zh) = 0; 

3. The vector meson spectrum should satisfy the linear Regge trajectories at zero temperature and zero 
density [36]. 
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With the above conditions, the equations of motion (2.8-2.12) can be analytically solved as 


b' (z) = + 


At (z) = fi 




1 - e“n 


(2.13) 

(2.14) 


9{z} 

Viz) 




2c//^ 
(l - 


2 


y^e ^"^dy y^e 
n„y^^~^^dy j^^y^e-^^e^y^dy 


—3z^ge 


-2A 


A" + 3A'‘^ + 



zj Z\2g 


zJ ^ 6g\ ’ 



(2.15) 

(2.16) 


where g = At (0) is dehned as chemical potential. 

The solution (2.13-2.16) depends on the warped factor A [z). The choice of A [z) is arbitrary provided 
it satishes the boundary conditions. To be concrete, we fix the warped factor A (z) in our solution in a 
simple form as 

Aiz) =-^z"^-bz^, (2.17) 

where the parameters b and c will be determined later. 


2.2 Phase Structure of the Background 

With the background (2.6), one can calculate the Hawking-Bekenstein entropy 


and the Hawking temperature 


^3A{zh) 


z%e 

rp _ M 


dvr y^e ^^dy 


1 - 


2c/i^ y^e ^^dy — y^e ^^e^y" 


1 - 


(2.18) 


(2.19) 


The temperature T v.s. horizon zh at different chemical potentials is plotted in figure 2. At /i = 0, 
the temperature has a global minimum T min (0) at zh = Zmm (0). The black hole solution is only thermo¬ 
dynamically stable for zh < -Zmin (0) and is unstable for zh > z^m (0). Below the temperature Tmin (0), 
there is no black hole solution and we expect a Hawking-Page phase transition happens at a temperature 
Thp (0) > Tmin (0) where the black hole dissolves to a thermal gas background. For 0 < y < the 
temperature has a local minimum/maximum Tmin (fJ-) /Tmax (l^) at zh = -^min (y) !Zmax id) and decreases 
to zero at a finite size of horizon. The black holes between Zmin id) and Zmax id) are thermodynamically 
unstable. There are two sections that are stable with zh < ^min (y) and zh > Zraaxid)- We expect a 
similar Hawking-Page phase transition happens at a temperature Thp id) ^ ^min id)- Nevertheless, since 
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T(GeV) T(GcV) 



Figure 2. The temperature v.s. horizon at different chemical potentials ^ = 0, 0.5, 0.714, IGeF. We enlarge a 
rectangle region in (a) into (b) to see the detailed structure. For /i > /ic, the temperature decreases monotonously 
to zero; while for /i < the temperature has a local minimum. At /ic — 0.714GeF, the local minimum reduces to 
a inflection point. 

the thermodynamically stable black hole solutions exist even when the temperature below Tmin for the 
section zh > f/r), we also expect a black hole to black hole phase transition happening at a tem¬ 

perature Tbb ifJ-) between Tmin (/r) and Tmax (fJ-), where a large black hole with the horizon z = zhi (fJ-) 
collapses to a small black hole with the horizon z = zhs (a*) as showed in figure 3. Finally, for /r > 
the temperature monotonously decreases to zero and there is no black hole to black hole phase transition 
anymore^. 

- z = 0 



Figure 3. Phase transition from a large black hole with the horizon z = zhi collapses to a small black hole with the 
horizon z = zhs at the transition temperature T = Tbb- 

To determine the phase transition temperatures Thp{^) and Tbs (/r), we compute the free energy 

^There could still be a Hawking-Page phase transition at some temperature for the case of /r > /ic, but we will show later 
that the black hole solution is always thermodynamically favored in this case. 
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from the first law of thermodynamics in grand canonical ensemble 



We plot the free energy v.s. temperature in (a) of figure 4. 



( a ) 


(b) 


( 2 . 20 ) 


J^^' ^(QeV) 
1.5 


Figure 4. (a) The free energy v.s. temperature at chemical potentials /r = 0,0.5,0.714, IGeF. At /r = 0, the free 
energy intersect with the x-axis at T = Thp where the black hole dissolves to thermal gas by Hawking-Page phase 
transition. For 0 < /i < /ic — 0.714GeF, the temperature reaches its maximum value where the free energy turns 
back to intersect with itself at T = Tbb where the black hole to black hole transition happens. For jj, > fj,c, the 
swallow-tailed shape disappears and there is no phase transition in the background, (b) The phase diagram in T and 
/i plane. At small /r, the system undergoes a first order phase transition at finite T. The first order phase transition 
stops at the critical point (/ic,Tc) = (0.714GeF, 0.528GeF), where the phase transition becomes second order. For 
> fic, the system weaken to a sharp but smooth crossover [20]. 

At /i = 0, the free energy intersect the x-axis at T = Thp (0) where the Hawking-Page phase transition 
happens. The black hole dissolves to thermal gas which is thermodynamically stable for T < Tpp (0). We 
fix the parameter b ~ 0.273GeP^ in Eq. (2.17) by fitting the Hawking-Page phase transition temperature 
Thp (0) with the lattice QCD simulation of Thp — 0.6GeV in [37]. 

For 0 < < He the free energy behaves as the expected swallow-tailed shape. The temperature 

reaches its maximum where the free energy turns back and intersects with itself at T = Tbb (z^) where the 
large black hole transits to the small black hole. Since the free energies of the stable black holes are always 
less than that of the thermal gas {Fgas = 0), the thermodynamic system will always favor the small black 
hole background other than the thermal gas background. When we increase the chemical potential h from 
zero to He the loop of the swallow-tailed shape shrinks to disappear at h = Z^c- For h > Z^c the curve of 
the free energy increases smoothly from higher temperature to lower temperature. 

The phase diagram of the background is plotted in (b) of figure 4. At // = 0, the system undergoes 
a black hole to thermal gas phase transition at T = Thp (0). For 0 < h < fJ-e the system undergoes a 
large black hole to small black hole phase transition at Tbb (z^)- The phase transition temperature Tbb (z^) 
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approaches to Tjjp at /r —>• 0 that makes the phase diagram continuous at /r = 0. The phase transition 
stops at n = He and reduces to a crossover for h > fJ-c- 

The phase diagram we obtained here in figure 4 is different from the conventional QCD phase diagram, 
in which crossover happens for small chemical potential and phase transition happens for large chemical 
potential. In [20], by comparing with the phase structure in lattice QCD simulation, the authors argued 
that this ’reversed’ phase diagram should be interpreted as confinement-deconhnement phase transition of 
heavy quarks in QCD. In this paper, we consider the same background as in [20] to study pure gluon QCD 
with one additional heavy flavour, and not light quarks. Since our model describes heavy quarks system in 
QCD, the flavour field in the matter action 2.3 should be associated to the mesons make up of heavy 
quarks, i.e. quarkonium states. By fitting the lowest two spectrum of quarkonium states mj/n = S.ODGGeD 
and = 3.685GeI4, we can fix c ~ l.lGGeV^ in Eq. (2.17). 

Nevertheless, there left a problem that, in the gravity side, it is commonly believed that the conhnement- 
deconfinement phase transition in the field theory side is dual to the Hawking-Page phase transition. 
Hawking-Page phase transition is the transition between black hole and thermal gas backgrounds. How¬ 
ever, in our gravity background, the phase transition is between two black holes for a non-zero chemical 
potential. Thus it is not consistent to consider a black hole to black hole phase transition in the gravity 
side to be dual to the confinement-deconfinement phase transition in QCD. In the following of this paper, 
by adding open strings in the background, we will study this issue more carefully to gain a more reasonable 
physical picture. 

3 Open Strings in the Background 

In this paper, we consider an open string in the above background with its two ends on the boundary of 
the space-time at z = 0. There are two configurations for an open string in the black hole background. 
One is the U-shape configuration with the open string reaching its maximum depth at 2 : = zq; the other 
is the straight configuration with the straight open string having its two ends attached to the boundary 
and the horizon at z = zh, respectively. The two configurations are showed in figure 5. Since the dual 
holographic QCD lives on the boundary, it is natural for us to interpret the two ends of the open string as a 
quark-antiquark pair. The U-shape configuration corresponds to the quark-antiquark pair being connected 
by a string and can be identified as a meson state. While the straight configuration corresponds to a free 
quark or antiquark. 

The Nambu-Goto action of an open string is 

Sng = j (3.1) 


where the induced metric 


Gab = 


( 3 . 2 ) 




x = 0 ^ ^ 2 






(b) 


Figure 5. Two open string configurations. In (a), a U-shape open string connects its two ends at the boundary at 
z = 0 and reaches its maximum depth at z = zq. In (b), two straight open strings with their ends connecting the 
boundary at z = 0 and the horizon at z = zh, respectively. 

on the 2-dimensional world-sheet that the string sweeps out as it moves with coordinates is the 

pullback of 5-d target space-time metric 

ds^ = —: 5 — ( g{z)dt^ -b ddP' -I—I > (3-3) 

2 V ayz) J 

where, to study the thermal properties of the system, we consider the Euclidean metric and identify the 
periodic of the time with the inverse of temperature as /3 = 1/T. 

3.1 Wilson Loop 

We consider a r x to rectangular Wilson loop C along the directions (t, x) on the boundary of the AdS 
space attached by a pair of the quark and antiquark separated by r. The quark and antiquark located 
at {z = 0,x = -^rl‘£) are connected by an open string, which reaches its maximum at {z = zq,x = 0 ) as in 
figure 6. 

It is known that taking the limit to /3 = l/T allows one to read off the energy of such a pair from 
the expectation value of the Wilson loop, 

(W(C)) ~ (3.4) 

where V (r,T) is the heavy quark-antiquark potential [26, 33]. 

In string/gauge duality, the expectation value of the Wilson loop is given by 

{W{C)) = I ~ p-Son-shell ^ (3 5 ) 
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Figure 6. Wilson loop as the boundary of string world-sheet. 


where Smg is the string world-sheet action bounded by a curve C at the boundary of AdS space and 
Son-shell is the on-shell string action, which is proportional to the area of the string world-sheet bounded 
by the Wilson loop C. 

Comparing with Eq. (3.5), the free energy of the meson is defined as 


V{r,T) = T Son-shell (r-,r). 

3.2 Configurations of Open Strings 

The string world-sheet action is defined by the Nambu-Goto action. 


(3.6) 


S = j (fiC = J d^^Vdit 


G, 


(3.7) 


where Gab = is the induced metric on string world-sheet. For the meson configuration, by 

choosing static gauge: = x, the induced metric in string frame becomes 


2A{z) 

taneb ^ 


= Gabdi°'d^° = —2 —9 i^) + 


,2A(z) 


v'2 


1 + 


9{z) 


dx^, 


(3.8) 


where the prime denotes a derivative with respect to x. The Lagrangian and Hamiltonian can be calculated 
as 


2A(z) 


C = VdetG = — 2 ~\/d (-2) + 


72 


9 {z) + z'"^ 


With boundary conditions 

z (^x = = 1 = 2 ) ~ = 0 ) = zq, z'{x = 0 ) = 0 , 

we obtain the conserved energy 


e2A(2o) _ 

n{x = 0 ) =- ^Vgizo)- 


(3.9) 

(3.10) 

(3.11) 

(3.12) 
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We can solve z' from Eq. (3.10) 


z' = 


where 


a{z) = 


{z) 

0-2 (^^o) 

(z) 


The distance r between the quark-antiquark pair can be calculated as, 


r = 


r2 /“^o 1 

J dx = 2 J dz— — J dz 


9{z) 


(zo) 


- 1 


(3.13) 

(3.14) 

(3.15) 


where zq is the maximum depth that the string can reach. The dependence of the distance r on zq at two 
different horizons are plotted in figure 7. 


r(fm) 



Figure 7. Separate distance r between quark and antiquark v.s. zq at fi = O.BGeV. 

We see that for a small black hole (large zj/), there exist a dynamical wall at Zm < zh where r' (zm) 
oo. The open string can not go beyond this dynamical wall, i.e. zq < Zm, even when the distance r 
between the quark-antiquark pair goes to infinity as showed in (a) of figure 8. While for a large black hole 
(small Zh), the open string can reach arbitrary close to the horizon, but there is a maximum value for the 
distance at r = xm- If the distance between the quark and antiquark is larger than vm, there is no stable 
U-shape solution for open strings so that the open string of U-shape will break to two straight open strings 
connecting the boundary at 2 : = 0 and the horizon at z = zh as showed in (b) of figure 8. 

In summary, for a small black hole, open strings are always in the U-shape; while for a large black 
hole, an open string is in the U-shape for short separate distance r < tm and is in the straight shape for 
long separate distance r > tm- Thus when a large black hole shrinks to a small one, we expect that a 
dynamical wall will appear when the black hole horizon equal to a critical value zh^ for each chemical 
potential /r, showed in figure 9. 
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Figure 8. (a) For a small black hole, open strings can not exceed the dynamical wall at z = Zm and are always in 
the U-shape. (b) For a large black hole, an open string will break to two straight strings if the distance between its 
two ends is larger than zm- 



Figure 9. When a large black hole with horizon zhi < zh^ shrinks to a small black hole with horizon zhs > zh/j., 
the dynamical wall at z = Zm appears when zh = zh^ ■ 

Since each horizon is associated to a temperatures for black holes, we define the transformation 
temperature corresponding to the critical black hole horizon at which the dynamical wall ap¬ 
pears/disappears, for each chemical potential /r. The dependence of on fi is plotted in figure 10. 

When there is dynamical wall, open strings are always in the U-shape. This means that quarks and 
antiquarks are always connected by an open string to form a bound state, i.e. the meson state in QCD. It 
is natural to interpret this case as the confinement phase in the dual holographic QCD. On the other hand, 
when the dynamical wall disappears, an U-shape open string could break up to two strings if the distance 
between its two ends is large enough. This means that a meson state could break to a pair of free quark 
and antiquark in QCD. We interpret this case as the deconfinement phase in the dual holographic QCD. 
Therefore, the transformation temperature is associated to the transformation between the confinement 
and the deconfinement phases in the dual holographic QCD^. 

^This transformation between confinement and deconfinement phases is not necessary a phase transition, it could be a 
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r(GeV) 



Figure 10. The temperature corresponding to zh^i., at which the dynamical wall appears/disaapears, for each 
chemical potential fj.. 


3.3 Heavy Quark Potential 

When the black hole horizon zh > zh^, there exists a dynamical wall at z = Zm < zh- Open strings are 
always in the U-shape and the heavy quark potential can be calculated as 


r2 r^o 

J dxC — J 


V = T = / dxC = 2 dz 


(j{z) 


1 - 


<7^ (^o) 

Cj2 {z) 


(3.16) 


-2 Jo \/g{z) 

In short separate distance limit r —)■ 0, i.e. zq 0, we expand the distance and the heavy quark potential 
at zq = 0, 


r-20 

r = 2 j dz 

'o 


•■L 

•J. 


V = 2 dz 


aiz) 


cr [z] 


- 1 


(^o) 

'l _ (^o) 

0 ““ L (^) 


= riZQ + 0{zl), 


V 1 

= — + 0 ( 1 ), 

Zo 


where 


n = 2 




du ( 


This gives the expected Coulomb potential at short separate distance, 


(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 


smooth crossover as we will show later. 
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where 


In long separate distance r —>■ oo, i.e. zq ^ Zm, we make a coordinate transformation z = zq — zqw^. The 
distance r and the heavy quark potential V become 


where 



fr (w) = 2zow 


/ 2 \ f - Zqw"^) 


fv (w) = 2 zqw 


1 - 


0-^ ( 2 : 0 ) 


a (zq - zon?) 

^/g {zo - zow"^) L" ( 2^0 - Zo'W^)_ 


(3.23) 

(3.24) 


(3.25) 

(3.26) 


We learn from figure 7 that the distance r is divergent at zq = Zm, and the same happens for the heavy 
quark potential. By carefully analysis, we find that this divergence is due to the integrands fr {w) and 
fr (w) are divergent near the lower limit u; = 0, i.e. z = zq ^ Zm- To study the behaviors of the distance 
and the heavy quark potential at zq = Zm, we expand fr {w) and fv (w) at tc = 0, 


fr (w) 

fv (w) 


2zo 


-2zog {zq) 


( y ' (20) ' 

(zo) _ 


+ 0{w), 


2zoa (zo) 


-2zog {zq) 


y ' (20) 
<y{zQ) 


+ 0{w). 


(3.27) 

(3.28) 


The integrals (3.23) and (3.24) can be approximated by only consider the leading terms of fr (w) and fr (w) 
near zq = Zm- This leads to 


r {zq) ~ Azq 


-2zog {zq) 


y ' (20) ' 
{zq) _ 


V (zo) cz 4 zo(T (zo) 


-2zog (zo) 


y' (zq) 
y{zo) 


a {zo)r {zq) . 


(3.29) 

(3.30) 


From the above expression, we obtain the expected linear potential V = agr at long distance with the 
string tension. 


dV 

dV / dzQ 

_ y' {zq) r {zq) + a {zq) r' {zq) 

dr 

. 0 =.^ d'^/dzQ 



(3.31) 


The temperature dependence of the string tension for various chemical potentials is plotted in figure 11. We 
see that the string tension decreases when the temperature increases. At the confinement-deconfinement 


- 14 - 




























Figure 11. String tension v.s. temperature at // = 0.5,0.678,0.714GeF. (a) The string tension decreases with 
temperature growing in the confinement phase, and suddenly drops to zero at T = The region closes to the 
transition temperatures is enlarged in (b) 


transformation temperature T^, the system transform to the deconfinement phase and the string tension 
suddenly drops to zero as we expected [38]. 

The behaviors of the heavy quark potential at short distance and long distance agrees with the form 
of the Cornell potential, 

V {r) =-—+ asT + C, (3.32) 

r 

which has been measured in great detail in lattice simulations 

Next, we would like to look at the r dependence of the heavy quark potential by evaluating the integral 
in Eq. (3.16), which is divergent due to its integrand blows up at z = 0. We simply regularize the integral 
by subtracting the divergent part of the integrand. 


Vr 



a{z) 


_V9 (z) 

[ (Y \ 


^[l + 2A'(0)z] , 


(3.33) 


where 

C (zo) = -h 4^4'(0) In zq. (3.34) 

^0 

After the regularization, we are able to calculate the heavy quark potential. The result is plotted in (a) 
of figure 12. For short separate distance, the potential is proportional to 1/r as expected. While for long 
separate distance, there exists a critical horizon zh/j, for each chemical potential /i. For small black hole 
with zh > ZHfj,, the potential is linear to r for r —>■ oo. While for large black hole with zh < zh^, the 
potential ceases at a maximum distance tm- Beyond rjvf, the U-shape open string will break to two straight 
shape open strings. 

It is helpful to use a sketch to describe the heavy quark potential and the phases transformation with 
temperature changing. We plot the heavy quark potentials at various temperatures at a fixed chemical 


- 15 - 
























(a) (b) 

Figure 12. (a) V v.s. r at fj, = Q.bGeV for zr > ZHfj., zh = zr^ and zr < ZRf^. (b) The sketch diagram for heavy 
quark potentials at various temperatures at a fixed chemical potential n > fic- 


potential^ /r > //c in (b) of figure 12. For T < T^, i.e. zr > zufi, the heavy quark potential is linear at 
large separate distance r with the slopes decrease with the temperature increasing. The linear potential 
implies that the system is in the confinement phase. While for T > T^, i.e. zr < zr^j,, the heavy quark 
potential admits a maximum separate distance rjvf, beyond which the open string breaks up to two straight 
strings and the total energy of strings becomes constant. The constant potential implies that the system 
is in the deconfinement phase. The confinement-deconfinement phase transform happens at T = T^. 


4 Phase Diagram 

In the previous sections, we have studied the thermodynamics of the black hole background. We obtained 
two black hole phases and studied the phase transition between them. We also added probe open strings 
in the background and studied their configurations of U-shape and straight shape, which correspond to 
confinement and deconfinement phases in the dual holographic QCD. In [20], we interpreted the black 
hole to black hole phase transition in the background as the confinement-deconfinement phase transition 
in the dual holographic QCD, leaving a puzzle that a black hole background does not correspond to the 
confinement phase in QCD in the original AdS/QCD correspondence. In this paper, we argued that the 
U-shape and straight shape of open strings should correspond to confinement and deconfinement phases in 
QCD, but the transformation between the two phases seems always smooth without phase transition. In 
this section, by combining these two phenomena, we are ready to discuss the full phase structure for the 
system of the open strings in the black hole background, corresponding to the confinement-deconfinement 
phase diagram in the dual holographic QCD. 

®For fi < ^c, the picture is similar but more complicated due to the black holes phase transition in the background. Here 
we just illustrate the general properties for the heavy quark potential and leave the details of the phase transition to the next 
section. 
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4.1 Confinement-deconfinement Phase Diagram 

Let us consider the configurations of the probe open strings first. We have found that for a small black 
hole with zh > zh^, open strings can not exceed a dynamical wall at z = Zm even the distance r between 
the quark-antiquark pair goes to infinity. This means that both ends of the open strings have to touch the 
boundary at z = 0, and the quark-antiquark pair is always connected by an open string in the U-shape to 
form a bound state, which corresponds to a meson state in the dual holographic QCD, as in (a) of figure 
8. We interpret this phase as the confinement phase in QCD. On the other hand, for a large black hole 
with Zh < ZHfi, the two ends of the open strings could also contact the horizon instead of the boundary. 
If the distance r between the quark-antiquark pair is large enough with r > tm, an open string of U-shape 
would break to two straight open strings as showed in (b) of figure 8. Thus the meson state would decay 
to a pair of free quark and antiquark. We interpret this phase as the deconfinement phase in QCD. We 
should remark that for a small black hole with zh > zh^i, it is impossible for an open string to break up 
due to the dynamical wall at z = Zm- Thus even in the black hole background, the holographic QCD could 
still be in the confinement phase. This clarifies the puzzle in [20] that a black hole background does not 
correspond to the confinement phase in QCD in the original AdS/QCD correspondence. The black hole 
phases, open string configurations and QCD phases are summarized in Table 1. 


Black hole 

String configurations for r —>■ oo 

Phase in QCD 

Small (zh > zhijl) 

U-shape 

Confinement 

Large {zh < zh^) 

Straight 

Deconfenment 


Table 1. Black hole phases, open string configurations and QCD phases. 

For each chemical potential /x, we have calculated the transformation temperature corresponding 
to the critical black hole horizon zh^- The result of v.s. // is plotted in figure 10. On the other hand, 
the phase transition temperature Tbb of black hole to black hole phase transition in the background was 
plotted in (b) of figure 4. To investigate the relationship between and Tbb, we plot both of them 
together in (a) of figure 13. We see that the two lines are close to each other but not exactly the same. 
The two lines intersect at (/ic,Tc) = (0.678GeU, 0.536GeU), where we define as the critical point‘d. For 
H < He when the temperature increases from zero, the black hole grows with the temperature and a phase 
transition eventually happens at T = Tbb{h) < where a small black hole with horizon zhs > zh^i 
suddenly jumps to a large black hole with horizon zhi < zh/j, as showed in figure 9. In the dual QCD, this 
implies that the confinement phase transform to the deconfinement phase by a phase transition. While 
for H > when the temperature increases from zero, the black hole horizon increases gradually with the 
temperature and continuously passes zh^ at T = < Tbb{h)- It means that the confinement phase will 

smoothly transform to the deconfinement phase as a crossover. Putting everything together, we obtain the 

"^We have defined (/ic,Tc) = (0.714GeF, 0.528GeF) as the critical values of the background phase transition in section 2.2, 
here we redefined them as the true critical values of the confinement-deconfinement phase transition. 
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final phase diagram for the confinement-deconhnement phase transition in QCD, plotted in (b) of figure 13. 
For the chemical potential less than the critical point /r < ^.c, we have confinement-deconfinement phase 
transition. While for the large chemical potential > fic, the confinement-deconfinement phase transition 
reduces to a smooth crossover. The critical point is located at fic = 0.678GeF and Tc = 0.536GeF. This 
result is consistent with the conclusion from the lattice QCD for the heavy quarks [37]. 


TiGeV) 




0.4 0.6 0.8 1 1.2 


^' H(Ger) 


( a ) 


(b) 


Figure 13. (a) The phase diagrams from the pure black hole background (red line for Tbb) and from the configu¬ 
rations of open strings (blue line for T^). The two lines intersect at the critical point locates at /ic = 0.678GeT and 
Tc = 0.536GeF. (b) The final confinement-deconfinement phase diagram. For /i < there is a phase transition 

between confinement and deconfinement phases (red solid line); while for /i > /ic, the phase transition becomes a 
crossover (blue dashed line). 


4.2 Meson Melting in Hot Plasma 

In the confinement phase, open strings are always connected in the configuration of U-shape, so that the 
quark-antiquark pair always form a bound state, i.e. a meson state in QCD. In the deconfinement phase, an 
open string could be in the configuration of either U-shape or straight-shape. For short separate distance 
f < tm, the open string is still in the U-shape. While for the long separate distance r > tm, the energy of 
the two straight strings is less than the free energy of an open string in the U-shape and the U-shape open 
string will break up to two straight strings, corresponding to that a meson state melts to a pair of free 
quark and antiquark. The phenomenon of mesons melting has been previously studied in [38-41]. In this 
work, we define the screening length as the maximum length r/vf, achieved by a pair of quark and antiquark 
in the bound state at a temperature T > Tbb{t) foi' T ^ Tc ot T > for /i > ^c- The screening length 
at a fixed chemical potential and temperature can be determined by the equation V — 2Fq = 0, where V is 
the heavy quark potential energy defined in (3.16) and Fg is the free energy of a straight string dehned as 

rzH p2A 

Jo ^ 
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We plot the ’melting lines’ of screening length versus chemical potential in figure 14. The screening length 
is a possible signal form Quark Gluon Plasma (QGP). Right after the collision QGP is formed and the 
temperature is high enough to be in the deconfinement phase. As temperature decreases, heavy quarks 
form bound states at melting temperatures higher than the deconfinement temperature. This means heavy 
quark bound sates can coexist with QGP. 




Figure 14. The ’melting lines’ for various screening lengths are plotted in (a) with tm = 
0.13,0.15,0.18,0.20,0.22,0.24 fm from above. In (b), the region closed to the phase transition line is enlarged. 


5 Conclusion 

In this paper, we considered an Einstein-Maxwell-scalar system. We solved the equations of motion to 
obtain a family of black hole solutions by potential reconstruction method. We studied the thermodynam¬ 
ical properties of the black hole background and found black hole to black hole phase transitions at the 
temperature Tbb for the backgrounds. We then added open strings in these backgrounds and identified 
the two ends of an open string as a quark and antiquark pair in the dual holographic QGD. By solving 
the equations of motion of the open strings, we got two configurations for the open strings, i.e. U-shape 
and straight-shape. When the temperature is low enough, the black hole is small, there exists a dynamical 
wall at z = Zm which the open strings can not exceed even the separation of the quark and antiquark 
goes to infinite. From the view of the dual QGD, the quark and antiquark pair is always connected by 
an open string to form a bounded state, corresponding to the confinement phase in QGD. On the other 
hand, when the temperature is high enough, the black hole becomes large so that an open string could 
break up to two straight open strings connecting the boundary and the black hole horizon, corresponding 
to the deconhnement phase in QGD. We obtained the conhnement-deconhnement phase transformation 
temperature T^. 

Our main conclusion is that, to study the conhnement-deconhnement phase structure in holographic 
QGD models, we need to combine two phase phenomena in the bulk gravity theory at the same time, 
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namely the black hole to black hole phase transition in the background and the various configurations for 
the probe open strings. We found that, when the chemical potential is less than the critical value //ci the 
background undergoes a small black hole to a large black hole phase transition with temperature increasing 
from zero. The horizon suddenly blows up to pass the critical horizon zh^i so that the confinement phase 
transform to the deconfinement phase by a phase transition. While when the chemical potential is greater 
than the critical value //c, the black hole horizon grows gradually and continuously pass the critical horizon 
ZHii so that the confinement phase transform to the deconfinement phase by a smooth crossover. The final 
confinement-deconfinement phase diagram is showed in (b) of figure 13. 

We also studied meson melting in this paper. When the temperature is higher than the phase transition 
temperature, QCD is in the deconfinement phase. However, it is known that the meson could still be 
thermodynamically stable in the deconfinement phase if the separate distance between the quark and 
antiquark is short enough. Only when the separate distance is longer than the screening length the 
meson becomes unstable and break up to a pair of free quark and antiquark. We showed the ’melting lines’ 
for various separation distance in figure 14. We conclude that, with increasing temperature, the mesons 
of larger size will break up earlier and the mesons of smaller size will be more stable and break up latter. 
Inversely, when QGP is cooling down, the mesons of smaller size will reunite earlier than the ones of larger 
size. This will help us to understand the process of the QGP cooling down. 
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